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We examine special extended spin 5 = 1/2 fermionic and hard core bosonic t — J z models with nearest 
neighbor and next nearest neighbor interactions to find exact ground states. Some of these models display an 
exponentially large degeneracy with diagonal stripe like patterns. 
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I. INTRODUCTION 

It has been nearly two decades since the discovery of 
the high-T c superconductors [j]]. Although much has been 
learned since, there is still no satisfactory explanation of what 
causes the superconductivity. There is widespread agreement 
that it should be possible to describe the electronic properties 
of the square CuC>2 lattices by the two-dimensional repulsive 
Hubbard model. This well known model is given by 



H 



where c\ a creates an electron on site i with spin a and j is a 
nearest neighbor of i. This model contains both the movement 
of the electrons (hopping) (i, kinetic energy) and the interac- 
tions of the electrons if they are on the same site (U, potential 
energy). The Hubbard model is one of the simplest possible 
models of interacting electrons. In the large U limit, the model 
reduces to the well known t — J model 



-* 2J 



(1) 



H 



E 



h.c. 



(ij) 
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and the now constrained fermion operators c| ,Ci ff obey 
the constraint of no double occupancy at any site i (rn = 



has expectation values or 1). In formal terms, 



the Fock space is projected to this sector adhering to these 
constraints. An even simpler version is afforded by the so- 
called t — J z model given by 



H = —t 



(ij),<- 



J 2 



2^ S < S 3 

(ij) 
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wherein the SU(2) spin symmetry of the original Hubbard 
and t — J models has been lifted. 

Numerical calculations on such models with open and 
cylindrical boundary conditions have found non-uniform 
"stripe" patterns [2] which have earlier been predicted by ap- 
proximate solutions to the Hubbard model | 3], and derived, 
very convincingly, by the addition of the strong Coulomb ef- 
fects present in the cuprates |4|]. At the moment, there is no 
clear consensus as to whether these patterns display the gen- 
uine ground state of the bare Hubbard (or t — J) Hamiltonian 



or amount to a finite-size artifact. Attacking these models 
has proven not to be an easy task. The Hilbert space of the 
model is too large for exact numerical solutions. Monte Carlo 
simulations suffer from the sign problem. The interaction en- 
ergy U is approximately equal to 8t, so we are neither in an 
extremely strong, nor in a weakly interacting limit. As has 
been emphasized by [4], Coulomb repulsions in the cuprates 
are poorly screened and universally promote the existence of 
stripes. Notwithstanding, it would be of interest to attain rig- 
orous results on whether various inhomogeneous electronic 
patterns can, in principle, be realized in the presence of short- 
range interactions and hopping alone in some models. 

The less ambitious goal of this article is not directly re- 
lated to the very challenging endeavor of understanding such 
complicated systems (albeit motivated by it). Rather, we con- 
struct exactly solvable high dimensional models related to the 
original t — J z models. Inspired by exactly solvable AKLT 
Uffl and in particular Klein models H, |1, H, 0, (H, 
we construct an exactly solvable model whose individual pair 
building blocks are an extension (containing superconducting 
pair terms and additional phases for the fermionic version) of 
the t — J z approximation of the Hubbard model. The price 
for attaining this solvable short range model (all interactions 
involve sites which are, at most, two lattice constants apart) 
and others is that their results cannot be directly related to the 
models of more physical interest (such as the plain t — J z 
model). 



II. CONSTRUCTION OF QUASI-EXACT SOLVABLE 
S=l/2 FERMIONIC AND HARD CORE BOSONIC 
SYSTEMS ON A SQUARE LATTICE 



In the up and coming, we follow a generalization of the 
Klein models H, H, Hi, O, (HI yet now for fermionic 
S = 1/2 particles on the square lattice. To do so, we first 
find, much in the spirit of Q], @|, 0, 0, [H a quasi- 
exactly solvable 5=1 model on a square lattice fl2il . We 
then apply a Jordan Wigner transformation to convert the S = 
1 problem into a fermionic 5=1/2 system. The resulting 
5 = 1/2 system is none other than a t — J z model augmented 
by superconducting pairing terms and phase factors. 
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A. Construction of ground states for the 5 = 1 system 

First, we generalize the Klein model (H, H, d, O, 
fiHo constructions to the problem of 5 = 1 spins on the 
square lattice whose sites are labeled by i (the locations are 
* = (ix,iy) are vectors with i x and i y integers). We assume 
that the lattice has periodic boundary conditions. For clarity, 
we dispense with the vector sign henceforth. We consider a 
spin 5 = 1 system whose Hamiltonian can be written as a 
sum over all plaquettes of the lattice 

H = Y J hn, (4) 

□ 

with ha a polynomial of the total spin of the four sites forming 
a given plaquette 

i.e. 

m 

= P(5^) = 5>„5* n n . (6) 

n=0 

Numerous polynomials (and hence short-range Hamiltonians) 
can be engineered to find exactly solvable models. Before 
doing so let us review the possible four and two particle states 
of 5 — 1 spins. With four spins in tow, the total spin of any 
plaquette, < 5p < 4. The Hilbert space spanned by four 
5 = 1 spins can be written as 

l®l(g)l<g)l = 

4® (3) 3 e(2) 6 ® (i) 6 e(o) 3 . (7) 

Here, the left hand side denote the direct product of the four 
5 = 1 spins while the right hand side denotes the decomposi- 
tion in the total spin basis. 

The sum of any two nearest neighbor spins on the lattice 
(Spair{S P ai7- + 1) = (5, + Sj) 2 ) can only attain the three 
values Sp a i r = 0, 1, 2. If any two spins forming a given pla- 
quette are in a singlet state {Spair = 0), then the total spin of 
the plaquette cannot exceed two, Sa _ 2. It follows that if the 
polynomial V of Eq.0 is non-negative definite and equal to 
zero for Sa — 0,1,2, then any state \ip) which is a superpo- 
sition of dimer coverings of the lattice each of which has (at 
least) one dimer per each plaquette then 

\i>)=Y J <*p n i<%>> w 

p ijeP 

(with such dimer coverings labeled by P) is a ground state of 
H. The above singlet state is given by 

\Sij) = -^\\ m * = = -!) 

+ | TOi = -1, m,j = 1) - | mj = 0, mj = 0)]. (9) 

Putting all of the pieces together, any state of the form 
Eqs.( l8l9> is a ground state of the Hamiltonian of Eqs. (0]|6j 



for any non-negative definite polynomial V which vanishes 
for Sa = 0, 1, 2. For concreteness, we now explicitly give an 
example of such a polynomial, 

ha = Sl ] (Sa-2)(S D -6). (10) 

Many other non-negative definite polynomials over all possi- 
ble plaquette spin values can be written down such that they 
vanish only for Sp = 0, 1, 2. Any of the dimer coverings of 
the form given by Eqs.( l8l9> is a ground state of the Hamilto- 
nian given by Eas.(l4|[T0t. 

B. Geometry of the ground states of the S = 1 system 

Enumerating all possible ground states of the form of Eq.(|8) 
amounts to the task of finding all possible dimer coverings 
of the on a square lattice with periodic boundary conditions 
such that every plaquette hosts, at least, one dimer. The prob- 
lem of hard core dimer coverings was investigated by Jl3ll 
wherein the number of fully packed coverings on an infi- 
nite lattice was found to scale as N ~ exp[GA/7r] with A 
the number of sites in the lattice and G Catalan's constant 
( G = E^o (L+i)* ~ 0.915965). When the condition of (at 
least) a single dimer per plaquette is imposed on a square lat- 
tice with periodic boundary conditions, we find that there are 
2 L (with L the lattice length) collinear diagonal coverings^. 
Along each diagonal, we may decide at our whim, to orient 
the dimers either vertically or horizontally. As the number of 
diagonals is L and as there are two such choices for dimer ori- 
entations (horizontal/vertical) for each diagonal, the number 
of individual dimer ground states scales as 2 L . Apart from 
these 2 L "stripe" states, there are no other dimer coverings 
such that one dimer appears in each plaquette. 

Implementing these geometric constraints on the general 
ground states of the 5 = 1 system of Eqs.( l8l9> , the allowed 
ground states read 

w=E«^in^(«.dbi.)>- 

P± ieA 

The difference between the general dimer states of Eq.(|8} and 
Eg. (II 1> lies in the explicit form of the allowed partitions P 
which is discussed below. In Ea. il 1> . the sites i £ A belong 
to one sub-lattice (say the one given by (i x ± i y ) being even 
with i XiV the x and y components of i) and the unit vectors 
4>(i x ± iy) = e x or — e y are random Ising variables. The 
configurations P + are those in which the orientation of the 
dimers (horizontal or vertical) is fixed along right-tilting di- 
agonals (i x — i y = const) whereas in P_ are those in which 
the orientation of the dimers is fixed along left-tilting diago- 
nals (fixed i x + i y ). As there are 2 L possible sequences in 
P+ for the field 4>(i x ± i y ) along all of the right-tilting diag- 
onals (whose number is L) with a similar number of possible 
coverings P_, the set of coverings {P+, P_} spans C(2 L+1 ) 
configurations. Similarly, additional dimer coverings where 
the left/bottom most position of the dimer may span both sub- 
lattices appear. Various orders can be stabilized by augment- 
ing the Hamiltonian by an additional field favoring this or the 
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other ground states |9[. A cartoon of such an external field 
of relative strength g stabilizing the parallel (all tj> = e x in 
Eci.i ll 11 1) and zig-zag phase (<fi = e x , e y , e x , e y , ... uniformly 
on consecutive parallel diagonals) at different values of the 
external field g is shown in Figure 1. Such a term may be, 
e.g., 



Si 



E $ 

ix+i y =o( mod 4) 

E 

n +i H =2( mod 4) 

-(S i + S i+ix ) 2 }]. 



S 



i+e y 



(12) 



Here, the first term of Eq. (ll2> . spans sites i belonging to half 
of the sites of sub-lattice A (the sites in which the sum of the x 
and y components of the sites i is an integer multiple of four, 
i x + i y = Am with m an integer). The second term in Eu. (ll2> . 
spans sites which belong to the other half of the sub-lattice A. 
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FIG. 1: Reproduced from @]- two possible dimer coverings of the 
square lattice (a, b) such that each plaquette hosts one dimer. In the 
context of the S = 1 Hamiltonian discussed here, the ovals denote 
a singlet state \Sij) of Eq.j5J formed between two nearest neigh- 
bor sites i and j. Such coverings can be stabilized by augmenting 
the Hamiltonian of Eq. J4l6> with additional terms similar to those of 
Eq. <12> . As a function of g a transition occurs between the single 
diagonal (SD) configuration of (a) and the zigzag (ZD) phase of (b). 
The shorthand O.P. denotes an order parameter in each of the phases. 

This model hosts a form of topological order in some of its 
ground states which is seen by string correlators. Any ground 
state 



n 

ijeP 



(13) 



is a direct product of two leg ladders (which are oriented along 
the diagonals) which are fully occupied by singlet states along 



each rung. Correlators involving spins along a line of consec- 
utive singlets in P alone display the same behavior as that of 
a two leg ladder. As the wave-function decouples into a prod- 
uct of two spin singlet wave-functions along each rung of the 
ladder, it is readily verified that in the states of Ea. (ll3t . the 
string correlator 



(Sf^exppTr ]T (Sf fc + S&fiSZjSZJ = -, (14) 



k=i+l 



irrespective of the separation between rungs i and j. Here, the 
indices 1 and 2 denote the location of the spin on one of the 
endpoints of the rung. A correlator of the form of Ea.(ll4> was 
first introduced by Todo et al. l H^ " 



C. An extended t — J z 5 = 1/2 fermionic model 

We now map our solvable 5 = 1 Hamiltonians and their 
ground states onto an electronic (5 = 1/2 fermion) problem. 
To this end, we employ a Jordan- Wigner transformation map- 
ping the three states of a 5 = 1 spin onto the three states of 
a constrained 5 = 1/2 fermion (one in which a forbidden 
double occupancy at any given site is strictly enforced). 

Rather explicitly, following 1 15], the mapping between 
these constrained 5=1/2 fermions and the 5 = 1 spin oper- 
ators is 



57 = V2(K 



5| = c 



:t 



c jt c J'T ~ c jl c Jl (15 ) 
Here the two dimensional kink operators leading to the correct 
(fermionic) statistics upon transmutation via Aharonov-Bohm 
phases by flux attachment 1 1 611 are 



K, 



exp[i 



E< 



k,jn k \ 



(16) 



with nk = cLcfc-r + cL Ckl the total occupancy of site k and 
9kj the angle between the line connecting sites k and j and, 
say, the horizontal x-axis. Here the sum spans all lattice sites 
k. 

Explicitly inserting this into the scalar products forming ha 
m Eqs.gl [10), we find that in ha any former appearance of 
(5 Q • Sb) in the 5=1 problem is now replaced by a t — J z like 
pair term augmented by superconducting pieces with phases 
in tow, 



t -,,_ e i4>ab j_-t £ i<t> a 



T a6 = 4^ + (cT T c 6T e^+cT T 



h.c. 



(17) 



Here, s ai b denote the 5 = 1/2 spin operators of the con- 
strained fermions on sites a and b which belong to the same 
plaquette, and the phases 



^6 



Ei 



EW 



hc)n c ] 



(18) 
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with the sum performed over all lattice sites c. The opera- 
tor T ab is an extended variant of the nearest neighbor interac- 
tion appearing in Eq.© with J z = At having additional terms 
(both additive and multiplicative). 

With this substitution any appearance of 5^ in the S = 1 
problem is simply replaced by 

S^2(^T a6 + 4) (19) 

(ab) 

where (ab) span all six pairs of sites which may be chosen 
from the four sites comprising the plaquette. The Hamiltonian 
governing the S —1/2 fermionic system is 

□ (ab)GD 

x( J2 T cd + 3)( Yl T f9 + 4)}- (20) 

(cd)en (fg)ea 

This is a particular variant of an extended t — J z model with 
fixed couplings. 

Similarly, the electronic state ground state of Ea. J20> (de- 
rived from \Sij)) is 

\xij) = A=Q Uj) + 1 UW - \0i0j)). (2i) 

Here, | denotes a single electron at site i with spin 

up (down) while |0j) corresponds to an empty site at site i. 
Thus a ground state of a Hamiltonian given by Eq.(I20> with 
the definition of Eq.ll^ is 

\i>)=Y, ap II ( 22 ) 

p ijeP 

In \ip), at all lattice sites k, the occupancy (nfe) =2/3 and 
the phases 4>ab,4>ab vanish if we replace, in Eq.( ll8t . by 
(nfe). Eq. J22t is superposition of states with different particle 
numbers; this is a consequence of the superconducting pair 
creation/ annihilation terms in Ea. (ll7> . The number of dimer 
coverings (the dimer partitioned denoted by P an example of 
which is shown in Fig.Q) is exponential in the perimeter of 
the system 

All stated above is exact. In the following paragraph (and 
here alone) we depart from the rigorous results derived above 
and briefly speculate about possible extensions. As within 
the ground state of Ea. (l22> there is no kinetic motion of the 
electrons- we might speculate that the transmutation statistics 
might well become irrelevant (as no transmutations occur). 
If this is the case, \ip) of Eq. (l22t . will be an eigen-state of 
Ea. (l20t with all phases <f> ab , <p ab set to zero. 

D. Hard Core Bose S = 1/2 extended t — J z models 

If we replace the hard core S — 1/2 fermions by hard core 
particles obeying other statistics, we find a generalization of 
the above. Rather explicitly, if hard core particles given by 



b a ,a,b aa , with a =f, 4, obey the fractional statistics of a 
phase change of e l7ra under the interchange of two particles 
(e.g. a = is the hard core bosonic case) then Ea. (l22> is a 
ground state of Eq.d20l subject to the substitution T a j, — * T" fc 
with 

T: b ee is* a st + (cl L c b] e^ +0^4/^ 

+c ai c fcT e- M *- + Za^er^ + h.c). (23) 

In the case of spin S — 1/2 hard core bosons [by which we 
refer to two flavor {a =|, J,) hard core bosons], the resulting 
Hamiltonian explicitly reads 

□ (ab)ea 

+c ai c 6T +c oJ cj 4 + h.c.) + 1) 

x( 4s cSd + (^jCdT+4r4i 

(cd)eD 

+ccicd<j + c c |4| + h.c.) + 3) 

x ( E 4s K + (4&i +z nii 

(fg)ea 

+c fl c g ^+c fl cl l + h.c.)+4,)]. (24) 

The states given by Eq. J22> are exact ground states of this 
short-range Hamiltonian. The frustrated kinetic hopping pa- 
rameters and competing nearest neighbor and next nearest 
neighbor antiferromagnetic interactions physically adhere to 
the resulting ground state. The precise form of the diago- 
nal formation of singlet dimers and high degeneracy of these 
states are non-trivial. Similar degeneracies appear in many 
other systems lfl8ll where effective dimensional reduction oc- 
curs. By extending the correlator of Ea. J14> to the 5 = 1/2 
arena by means of the transformations of IU5II and noting that 
exp[27ris|] = exp[i7m.;] at all lattice sites i, we find a coun- 
terpart to the Ogata-Shiba correlators fPUl for two leg ladders 

(xp\4,i4,i 

exp[i?r Y ( n M + n 2,k)}sf,j s 2, 3 \xp) = gg- (25) 

k=i+l 

Here, \xp) = U lj£ p \Xij)- 

Other ground states can easily be proven along the same 
lines. As the fully ferromagnetic states Sf — +1 for all i or 
Sf = —1 are ground states of the spin 5=1 Hamiltonian 

□ 

the maximally occupied (by hard core bosons) fully polarized 
ferromagnetic states sf — 1/2 (rij = 1) for alH or s\ — —1/2 
(jii = 1) at all sites i are ground states of the s = 1/2 Bose 
extended t — J z model 

H = - X)[ J2 4s >b + + clAl 

□ (afc)GD 

+Caic 6T + c al c f bl + h.c.)}. (27) 
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FIG. 2: (Color online.) From HHl . Highly regular ground states 
on the checkerboard and pyrochlore lattices. The ovals denote sin- 
glet dimer states. The arrows denote the representations of these 
dimer states within the six-vertex model (see text). On each plaque- 
tte (tetrahedron) the dimer connects the bases of the two incoming 
arrows. 



m. OTHER LATTICES 

We may repeat our derivations above to derive the ground 
states on any other lattice in which four fermions in any given 
basic unit interact with each other according to Eq.j20L One 
such example is afforded by the three dimensional pyrochlore 
lattice and its dimensional checkerboard rendition. The py- 
rochlore lattice is composed of corner sharing tetrahedra such 
that each vertex is common to two 4-site units (see the right 
hand side of Fig.©. The geometry is rather common in tran- 
sition metal and rare earth oxide systems, specifically in the 
pyrochlore (the origin of the name) and spinel structures. The 
checkerboard lattice (see the left hand side of Fig.©) is a 2d 
projection of the pyrochlore structure: the crossed plaquettes 
in the left hand side of Fig.© represent the basic tetrahedral 
units. 

On the pyrochlore/checkerboard lattices the open structure 
of tetrahedra/crossed plaquettes creates, by comparison to the 
square lattice, a less constrained system, and the number of 
dimer configurations in the ground-state manifold is strongly 



enhanced, the still more massive degeneracy implying simi- 
larly exotic physics in this case. 

In addition to the simple coverings shown in Fig. |2] a far 
richer variety of states exists. A detailed study of dimer states 
on the pyrochlore and checkerboard lattices is provided in 
Holl . A mapping onto the spin-ice problem provides a useful 
classification of these ground states. Here, the states with one 
dimer per tetrahedral unit defined by Eq. © map exactly to 
the spin-ice problem lfl9t l20ll where, motivated by the struc- 
ture of H atoms in solid water two of the sites of any 
elementary unit (a tetrahedron) are associated with an ingo- 
ing arrow pointing towards the center of the tetrahedron and 
two sites lie on arrows pointing outwards. We label the two 
sites belonging to a singlet dimer |SV, ) in a given tetrahedron 
□ by two incoming arrows from sites i and j to the center of 
the unit □ = ijkl. In this fashion it is clear that the system 
is mapped to a set of continuous directed lines such that each 
tetrahedral unit has exactly two incoming and two outgoing 
arrows relative to its center. An example of this mapping is il- 
lustrated in Fig. |2 The mapping is one-to-one: any spin-ice 
configuration determines a unique singlet-covering state with 
one dimer per tetrahedral unit and vice versa. The algebraic 
correlations in the spin-ice model dictate an algebraic decay of 
correlations in the d — 2, 3 dimensional 5 = 1/2 Fermi prob- 
lems of Eq. (I20> on the checkerboard and pyrochlore lattices : 
(T a bT mn ) ~ \r\~ d for large separations |r|. Here, |r| denotes 
the distance between the two site pairs ({ab) and (mn)) d, 



IV. CONCLUSION 

In conclusion, we constructed spin S = 1/2 extended t— J z 
models (both fermionic and hard core bosonic) which harbor 
easily provable exact ground states. We illustrated how this 
method can be used to engineer non-uniform diagonal dimer 
state and other configurations. Applying similar ideas, many 
other solvable systems can be constructed on other lattices. 

This work was partially sponsored by the Swedish Research 
Council. 
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